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Introduction 



>. 

CO . Fracture functions play a fundamental role in the description of one particle inclusive deep inelastic scattering 
■ H E H 0. They have been incorporated into the analysis of a variety of specific processes, although only up to 
\ leading order accuracy |E IE 0, since their scale dependence, which includes non homogeneous effects, had 

not been studied at higher orders until recently. At the lowest order in QCD perturbation (LO), these functions 
account for processes in which the final state particle is produced in the direction of the target remnant, and can be 



interpreted as the probability to find a parton in an already fragmented nuclcon target 0| . Their non homogeneous 
scale dependence is found to be strongly suppressed at LO in most of the kinematic domain, particularly in the 
kinematic region covered by the present experiments, and thus neglected. However, there are no indications if this 
suppression is also present at NLO level. 
CIh In reference |lO| the gluon initiated corrections to one particle inclusive deep inelastic scattering were computed 
^ up to order a1 for the first time. There it was shown that at variance with the fully inclusive case, in the one 
particle inclusive cross section at this order, coUinear configurations lead to singular distributions in more than one 
variable which can not be handled with the standard procedure. However a suitable method for the prescription of 
the characteristic entangled singularities was developed, verifying explicitly the factorization of coUinear singularities 
and obtaining the corresponding non homogeneous evolution kernels. Regarding the phenomenological consequences 
' of these corrections, it was found that the non homogeneous contributions to the evolution equations at NLO may 
be noticeably larger than the ones found at LO. Of course, gluon initiated corrections do not provide all the kernels 
relevant for the scale dependence of fracture functions, a consistent assessment of this dependence requires also the 
consideration of quark initiated processes at the same order. 

In the present paper we extend the computation of the ©(a^) corrections to one particle inclusive DIS for the case 
of quarks in the initial state. These last processes involve two caveats, a more complex singularity structure in the 
partonic cross sections due to the presence of virtual corrections, and complications associated to the flavor structure 
and decomposition for these corrections. 

Quark initiated corrections imply the computation of several real and virtual contributions which lead to a much 
more complicated singularity structure than in the gluon initiated case. The entanglement of singularities is precisely 
the most difficult point in the computation of higher order corrections to one particle inclusive processes. Nevertheless, 
we show that they can be handled extending the strategy implemented for the latter case. 

Regarding the fiavor structure, the explicit factorization of coUinear singularities in the one particle inclusive cross 
section does not allow to obtain all the evolution kernels individually, only certain combinations of them can be 



•Partially supported by CONICET, Fundacion Antorchas, UBACYT and ANPCyT, Argentina, 
t Electronic address: |daleo<!9tisica. unip . edu . ar | 



2 



extracted. Consequently, at this point it is not possible to achieve a complete flavor discrimination in the evolution of 
fracture functions. However it is possible to define certain singlet and non singlet flavor combinations, identify their 
respective evolution kernels, and obtain the evolution of the cross section. 

As in the case of gluon initiated corrections, NLO contributions to the non homogeneous evolution initiated by 
quarks are found to be comparable or even larger than the LO ones in some kinematical regions. In order to assess the 
significance of non homogeneous contributions relative to homogeneous ones in the scale dependence, we implement 
the full evolution equations for Mg .^+ /p using as input a model estimate for these fracture functions at an initial 
scale. The results we found indicate large non homogeneous effects at small xl and xb but also suggest suppression 
for both LO and NLO contributions otherwise. On the other side, this estimate does not take into account the NLO 
non homogeneous kernels for gluon fracture functions, which show up at order in the cross section. These kernels 
contribute indirectly to Mq through the homogeneous terms in the evolution but have not been computed yet. 

In the following section we summarize the relevant kinematics and details about the phase space integration for 
the 0{al) contributions to the cross section, together with the conventions and notation adopted. In section II we 
compute the corresponding real and virtual amplitudes associated to each of the quark initiated processes, we discuss 
the nature of the singularities that contribute to them and the prescriptions required. Then, in section III we analyse 
the factorization of coUinear singularities, and issues related to the flavor decomposition. There, we also present the 
evolution kernels obtained and derive the evolution equations. Section IV discuss the phenomenological implications 
of the new corrections and finally we end with our conclusions. 



I. NOTATION AND KINEMATICS 



We begin by reviewing the main features of the QCD description of one particle inclusive deep inelastic scattering 
processes. We consider the process 

l{l) + P{P)^l'{l') + h{Ph) + X, (1) 

where a lepton of momentum / scatters off a nucleon of momentum P with a lepton of momentum /' and a hadron h 
of momentum Ph tagged in the final state. The cross section for this process can be written as 

^ — ^ / HI 11 ■ I ^ ' \ II I ' \ n- I rinr T-> fi'ii fin : rill) ■ 



dxB dy dvhdwh ^ j^^^ gJxB ^ '^j Jo V ^ / / dxs dy dvj dwj 



E /' - (— - (1 - ^bH) (I-xb) - wn) , (2) 



where the sum is over all parton species. The leptonic final state, in the one photon exchange approximation, is 
described by the usual deep inelastic variables 

Q' = -q' = Hl'-ir, ^B^^, y^^, Sh^{P + 1)\ (3) 

2P ■ q P ■ I 

whereas, to describe the hadronic final state, we define 

Eh l-cos6h ... 

Eq(1-xb) 2 

Efi and Eq in the preceding formulae are the energies of the final state hadron and of the incoming proton in the 
P + q = frame, respectively. 6h is the angle between the momenta of the hadron and that of the virtual photon in 
the same frame. The convolution variables Vi and Wi in equation ^ are the partonic analogs of Vh and Wh and u is 
related to the fraction of momentum of the incoming parton ^ by ^ = xb/u. 

The two terms in the r.h.s. of equation (01 stand for 'current' and 'target' fragmentation processes respectively. 
In the first one d&ij represents the partonic cross section for the process l + i—^l'+j + X, whereas fi/p and D^/j 
are the usual partonic densities and fragmentation functions. In the second term the cross sections for the inclusive 
process I + i — > I' + X, d&i, are convoluted with the fracture functions Mi ^/p which give the probability of finding 
parton i in an already fragmented target. In this last case, the hadron is produced in the direction of the incoming 
proton (notice the (5(1 — Wh) in the last factor of the second term). 

The partonic cross sections d&i and duij can be computed order by order in perturbation theory from the corre- 
sponding parton-photon squared matrix elements, i + 7 — > X and i + 7 — s- j + AT respectively. The totally inclusive 
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cross sections up to order c?^ , are well known and can be found in reference ■ The lowest order one-particle inclusive 
cross section (zeroth-order in a^) is given by 

- - {2 + e)YM5{l-u)5{l-v)5{w), (5) 



dxB dy dv dw dxB dy dv dw 2 xb Sh 

where the outgoing quark is always produced backwards (w = 0). The index M stands for metric, longitudinal 
contributions are absent at tree level, and Ym = (1 + (1 — ?/)^)/2j/^. The result correspond to d = 4 + e space-time 
dimensions. 

The order-Qfj corrections to the one-particle inclusive cross sections are also known. Expressions for the singular and 
finite terms computed within the framework of dimensional regularizationj^^ can be found in 2] for the unpolarized 
scattering, and in Q for the polarized case. Finite contributions up to order e, needed for the factorization of collinear 
singularities at 0{a'^), are given in 

At order-Qfs the fragmenting parton can be produced in any direction, including the forward one {w = 1), but all 
singular contributions come either from the backward or from the forward direction. The former terms are factorized 
into partonic densities and fragmentation functions, whereas the forward singularities can only be factorized in the 
redefinition of fracture functions. This factorization gives rise, as we have already mentioned, to non-homogeneous 
terms in the evolution equations of fracture functions l|. 

At order-Q!^ the partonic cross sections develop singularities in every direction in space. They receive contributions 
from both virtual and real amplitudes whose matrix elements have to be integrated over the phase space of the 
unobserved partons (i.e. those partons that are not attached to the fragmentation function). At this order, real 
contributions have three particles in the final state, and the corresponding phase space can be written as: 



1 (Q'X (\_xb\^^^Uu-x_b 
(47r)4 r(l + e) Utt/ V a;s y ^ xb 

X (1 - wf^'^dv dw sini+" /3i sin' 02 dPi d/Ja . (6) 

The angles /3i and /32 are the polar and azimuthal angle of the pair of unobserved partons in their center of mass 
frame. In order to enforce the kinematical constraints and deal with the singularities in the prescription stage, it is 
convenient to break up the integration range into three regions, i? = BO U Bl U B2 with 

BO = {u G [xb,Xu], V G [vh,a\, w G [0, 1]} 
Bl = {u G [xb,Xu], V G [a, 1], w G [0, w^]} 

B2 = {u G [xu, 1], V G K, 1], w G [0, Wr]} , (7) 



where 



Xb _{1-u)xb _ {l~v)il- u)xb 



{xb + {I - XB)vh) ' u{1-xb)' ' v{u-xb) 

Bl and B2 receive contributions from both 0{as) and 0{a1) processes, while BO has only 0{ag). Virtual contributions 
can have either one or two particles in the final state. As it will be discussed in the next section, only diagrams with 
two particles in the final state give rise to contributions in the forward direction. In this case, the phase space can be 
conveniently cast into 

SttI (l-fe/2) \ 4:71 J u — XB \ xb J 

where now the integration is restricted to the region 1/ = Bl U B2. Notice that, for these contributions, the energy 
and angle variables, v and w, are correlated by the 5 function which constraints the variables to the surface w = Wr- 
This distinctive feature of the two particle phase space already appeared at 0{as) 01- 
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II. ORDER PARTONIC CROSS SECTIONS 



Order-a^ contributions to the one-particle inclusive cross sections come from the following partonic reactions: 



Real contributions 



Virtual contributions 



1 + 9 

J + q{q) - 

1 + qiq) - 

1 + 9 



> 9 + 9^ 

^ qt{qi) ^ 
^ qi{qi) ^ 

> 9 + q^ 

9i9) 

9 + q{q) 

q + q 



qj + qj 



qi 



(10) 



where any of the outgoing partons can fragment into the final state hadron h. Gluon initiated reactions were already 
discussed in deep in reference [lol |. In this section we analyze quark initiated processes, and examine the the nature of 
the singularities that they involve. These contributions are computed in d = 4 + e dimensions, in the Feynman gauge, 
and considering all the quarks as massless as in reference [Tni l . Algebraic manipulations were performed with the aid 
of the program Mathematica and the package Tracer to perform the traces over the Dirac indices. 

Taking into account all the alternatives for one parton undergoing fragmentation into the final state hadron h, the 
reactions in equation (|1U|I lead to the following partonic cross sections: cr^g, dq.q-, aq.g^, cTq.g-, cTq.g-, and their charge 
conjugates, exhausting all the possible combinations. The first index labels the parton that initiates the reaction, while 
the second stands for the one that precedes the final state hadron. Notice that for the last four cross sections, with 
quarks in both the initial and the hadronized state, the are several flavor combinations, each leading to characteristic 
singular contributions. 

Contributions to the aqg cross sections come from the first and the sixth reactions in equation (|10|l . corresponding 
to real and virtual diagrams respectively. These are shown in figures ^ and |2] Virtual contributions are obtained from 
the interference of diagrams in figure 12 with those of single bremsstrahlung at order as. The diagrams in the second 
row of Figure ^ need to be taken into account due to the symmetric character of the outgoing gluons wave function 
which also contributes with a factor 1/2 to the squared amplitude. 







FIG. 1: Real contributions at order for the 7 + q{q) g + g + q{q) process 



For the remaining contributions, where a final state quark hadronizes, it is necessary to discriminate whereas this 
last quark carries the same flavor of the initial state quark or not. Contributions to the <Tq^q^ and cTq-q^ cross sections 
with i ^ j, arise form the diagrams in the first row of Figure |3| corresponding to the second reaction in equation 
(|10|l . The amplitudes of the first two diagrams in this row are proportional to e^, whereas the two last amplitudes 
are proportional to Cj. Consequently, the above mentioned cross sections can be decomposed into three pieces, 
proportional to e| , e| and ej respectively. Since these cross sections first appear at order-a^ they receive no virtual 
contributions at this order. As the amplitudes are symmetric under the interchange of the quark-antiquark pair of 
fiavor J, the relation Oq-q. = cfq^qj is satisfied. 

Finally we have the contributions to (Jq^q^ and (Jqiqi^ where the initial and hadronizing quarks carry the same flavor. 
In this case, the relevant reactions are the first three, the fifth and the sixth in equation (|10|l corresponding to diagrams 
in Figures n and 131 (contributions in both i ^ j and i = j cases) for the real contributions and diagrams in Figures 
121 auditor the one and two loop virtual corrections respectively. Notice, that in spite of having the same label for 
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FIG. 2: Virtual contributions at order for the 7 + q{q) g + q{q) process 




FIG. 3: Real contributions at order for the 7 + q{q) q + q' + q{q) process. Diagrams in the second row only contribute 
in the case q = q' 

initial and hadronizing quarks, these cross sections include also terms proportional to ef, '^j^^^'j and '^j^^eiej, 
since the quark which enters the electromagnetic vertex in the third and fourth diagrams in figure |31 can be of any 
flavor. Regarding the process j + Qi ^ Qi + Qi + Qi, diagrams in the second row in Figure |31 have a relative minus sign 
respect to those in the first row due to the antisymmetry of the wave function of the identical outgoing quarks which 
also contributes with a factor 1/2 to the matrix element. 

The angular integration of the matrix elements for the real contributions can be performed with the standard 
techniques |l5l |. taking into account the additional complications of the one particle inclusive case: the necessity of 
collecting to all orders the potentially singular factors in the three particle final state integrals, as explained in 10]. 
For the integrals that are known to all orders in e, this is not a problem, while for those which are only known up to 
a given order a careful treatment is required. In Appendix A we include a list with the integrals that appeared in the 
present calculation. 

Once the angular integrals are performed, matrix elements are still distributions in u, v and w, regulated by the 
parameter e. In order to write these contributions, specifically their potentially singular pieces, in a consistent way, 
they need to 'prescribed', i.e. written as an expansion in powers of e with integrable coefficients. At variance with 
the totally inclusive case, where the integration over final states leads only to singularities in the u variable, in this 




FIG. 4: Two loops contributions at order for the 7 + q{q) q{q) process (quark form factor). 

case they appear along various curves in the residual phase space. Contributions with singularities along one or 
two intersecting curves can be managed with the technique developed in 10] . However, processes with identical 
quarks in both initial and hadronizing states lead to singularities along three intersecting curves, requiring a different 
prescription strategy. This distinctive feature of the aq-q- cross section is related to the fact that it give rise to quartic 
poles in e, which cancel out between real and virtual contributions. 

As an example, we consider the subprocess j + q—^q + g + g, where the final state quark fragments into hadron 
h. The quark propagators coming from the interference of the first and fourth diagrams in the first row of Figure ^ 
contribute with a factor 



1 



(fcl +p/)2 (p, - k2f (fcl + k2 +PfY {p^ - fcl - fc2)2 



(11) 



where pi and pf are the momenta of the incoming and outgoing quarks respectively, while fci and fc2 are those of the 
gluons. In the center of mass frame of the two gluons, the integration over their phase space involves only the first two 
denominators. In this frame, the two other denominators do not depend on the angles f3i and (32, defined in equation 
©. In this particular case, the angular integrals can be performed to all orders in e and have single poles in e = 0. 
The most singular terms can be written as 



1 

— w 

e 



l+e 



(W-Wo) ^' (Wr-W)' , 

1 — u 



(12) 



where /(u, v, w) is a regular function in all the integration region, specifically finite when e ^ and Wq = —(1 — v)/v. 
Notice the presence of singularities along w — 0, w — wq and u = I, and that these three curves intersect at the point 
u — V = 1, w = 0. 

In region BO only the factor w~^^'^ in equation H12I) is potentially divergent, since v < a < 1 and u < x„ < 1. This 
allows the use of the standard prescription recipe, which lead to a double pole accompanied by a S{'w). In region 
Bl, the singularity at w = wo need also to be taken into account, however this can be managed with the approach 
described in this case we obtain 



w) 



= /2 



Bl 



. 3e/2n we/2 1-2. r(l + e) r(l + e/2) 

-<5Ma/(l-a) r(l + 3e/2) 



e:2e, 1 + -e;a 



xi'^^-^^(l--) + ((l--)-^"")..[.^ 



2e 



-\-w [0,1^1-] 



(13) 



Making this substitution in equation (|12l) . we find that, in region Bl, the leading singularity is a triple pole, propor- 
tional to S{w) (5(1 — v). 

Finally, in region B2, the factor (1 — u)^^ is also potentially divergent, and should be taken into account. In order 
to do that, we use the prescription given in equation H13|l replacing a with z as the lower limit in the v prescription. 
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obtaining 

£ du 1^ dv j^^ dw ^5{w) (1 - u)-^+^ [^1Z^5{1 ^v) + ((1 - «)^^+'^)+,[,^^) ~f{u, V, w, e) 

+ [ du [ dv [ dw-^— (w-^+^iw-way^+'/^iwr-wy/A f(u,v,w,e), (14) 

Jx^ Jz Jo e 1-u \ /+wl0.w,.] 

In the last equation we have factored out in f{u, v, w, e) the expressions that are not relevant for the prescription. In 
the first term, the factor (1 — u)~^~^'^ can be prescribed with the usual recipe 

(1 - -^-jii-u) + ((1 - ur'+^)^^^„^^ , (15) 

whereas the second term can be safely expanded in power series of e, since the u integral is convergent even when 
e — > 0. Indeed, in this case, the pole in m = 1 is regulated by the fact that the w integration region shrinks to the point 
w = when u 1 whereas the integrand has no 6 functions. The final result includes a quartic pole, as mentioned 
above, accompanied by S{w)S{l — v)S{l — u). 

All the terms singular along three intersecting curves found in the quark initiated corrections can be managed as in 
the previous example, and they present the distinctive feature of having singularities along u = 1. Consequently, as 
it is shown above, the corresponding prescriptions have to be modified, but only in region B2. Given that the NLO 
non homogeneous kernels are determined by singularities along w = I exclusively, associated to the prescription in 
the BO region, and to terms which contribute in the forward direction in Bl, they receive no contributions from the 
more involved triple overlapping singularities. 

The virtual one loop contributions to both (jqg and (Tq^g. can be straightforwardly computed using the standard 
techniques for handling the loop integrations. These corrections, with two partons in the final state, contribute 
to regions Bl and B2 but not to region BO. In this case, the singularities relevant for the renormalization of fracture 
functions, that is in the w = 1 direction, are along the curve v = a in region Bl and can be managed with the usual 
prescription recipes in one variable. On the other side, due to the one-body phase space of the outgoing quark, the 
two loop graphs in Figure01only have support in the point u = v — l,w = and thus their calculation is not required 
for the extraction of the non homogeneous kernels. The results for these contributions can be found in [TtL lla. IT^. 
Of course, the complete calculation of the 0{a^) cross sections requires to take into account not only the two loop 
quark form factor but also a full analysis of the singularities coming from the real and one loop virtual contributions, 
including regions Bl and B2. 



III. FACTORIZATION OF SINGULARITIES AND EVOLUTION EQUATIONS 

The order-a^ one particle inclusive partonic cross sections described in the previous section, have a complex singu- 
larity structure, revealed as poles in e. Poles associated with ultraviolet divergences are removed by coupling constant 
renormalization. On the other side, factorization theorems warrant that poles related to coUinear singularities can 
be removed from the hadronic cross section by suitable redefinition of the bare parton densities, fragmentation and 
fracture functions in terms of renormalized ones. The self consistency of these factorization prescriptions fixes the 
scale dependence of the renormalized quantities, and in this case allows to obtain the NLO evolution kernels for the 
fracture functions. 

In order to accomplish the factorization procedure, we begin writing the hadronic cross sections in terms of renor- 
malized parton densities and fragmentation functions. Restricting our attention to the region BO, which contain all 
the forward singularities, we avoid the appearance of divergences at u = 1 or u = 1 as explained in Section II. Notice 
that in this way, we only have to take into account LO expressions for the bare parton densities and fragmentation 
functions, since NLO contributions, once convoluted with the tree level partonic cross section, give rise to terms with 
support in u = 1 and v — 1 which contribute only to regions Bl and B2. 

As explained in reference jlflj, the convolutions between the LO kernels and the 0{as) cross sections can be better 
handled keeping to all orders in e the expressions for the cross sections, and writing the subtracted 'plus' distributions 
in the kernels as 

(t^) - lim (1 - x)-'+^' -^6{1- x) + 0{e') . (16) 

The limit e' — > can be safely taken once the convolutions are explicitly evaluated. The resulting counterterms can 
then be subtracted from the 0{a'j.) cross sections while ultraviolet singularities are removed by coupling constant 



8 



renormalization 

e/2\ 



as _ a, 
2^~~ 



{Ml) as {Ml) ^ PofMlY 

27r I 27r ''^ e J 



(17) 



/3o is the lowest order coefficient function in the QCD P function, 



/3o = y - ^ UfTp (18) 

where Ca = N for SU{N), Tp = 1/2 as usual, and n/ stands for the number of active quark flavors. Mr is the 
renormalization scale and 

/r = ^^^^''^^ (19) 

r(i + e) ^ ' 

At this point, only poles proportional to 5{1 — w) are left in region BO, to be factorized in the redefinition of the bare 
fracture functions. In terms of renormalized quantities, fracture functions can be written as 

1 y A.^, (z., M,) m;:,/p (i C, M|) . (20) 

where a sum over repeated indices is understood. Notice that we have chosen all the factorization scales to be equal 
to the renormalization scale. 

Given that only single and double poles remain in the cross sections once that coupling constant, parton densities 
and fragmentations functions are renormalized, the functions A^^j and A^i^j can be written as 

f2 \ n , _ . 9 / A#2 ' 



The Ai,_j functions in the homogeneous terms are just those found for parton densities in the inclusive case, and 
can be obtained from the corresponding transition functions in reference ^llj . The LO non homogeneous kernels 
P^^^ (u,w), can be written in terms of those in reference |3] as 



P^, {u, v) = Pj^l^ (u)S[v~—^] . (23) 

Whereas the NLO kernels Pj^ll^j (u, v) can only be found after explicit calculation of the partonic cross sections 
up to order a^, the requirement of consistency in the factorization process constraints the double pole coefficients 
i?^^^^-(u, u). Indeed, in order to factorization be consistent, equation H20|) has to be valid no matter the election of 
the factorization scale M f and the corresponding renormalized quantities have to be well defined when e ^ 0. These 

conditions imply that the 0{a^) coefficients -R^^^j (u, w) are determined by the 0{as) kernels. Explicitly, 
where the convolutions are defined as 



1 -It 



f{u,v)®g{u) = / ^ f{u,v)g (^^^ , 
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f{u,v)^g{v) = / " , (25) 

f{u,v)(g)' g{u) = [ ^^f(u,^v)g(^). 

/,,, u u V u / \u/ 



1-"^ duu / u 
— - f [u, -1 

The fact that equation H24|l actuaUy holds for any combination of indices is a useful mean of checking the results and 
provides a non-trivial verification of factorization of coUinear singularities. 

Once the bare fracture functions in the hadronic cross section, equation are written in terms of the renormalized 
quantities, equation H20I) . the explicit cancellation of forward coUinear singularities, showing as single poles in the 
partonic cross sections, define the NLO evolution kernels for fracture functions. 

Forward singularities coming from the aqg cross section determine the non homogeneous evolution kernel 

Pgl}-q{u,v). Notice that the first index to the left in this notation refers to the parton that fragments into the final 
state hadron, the second labels the flavor of the fracture function that factorizes the corresponding pole, and the one to 
the right is for the initial parton. Due to charge conjugation symmetry aqg = a^g and thus PggLg(u, v) — Pgq]^q{u, v). 
The kernels coming from the crg.^^. cross sections can be classified identifying the cases in which the quark that 

interacts with the photon carries either the flavor of the initial state or that of the quark that hadronizes, Pq^l-^q. {u, v) 
and Pq^^.^q.{u,v) respectively. Notice that in the first case the corresponding pole in the cross section will be 
multiplied by ef , in the second by ej, while there are no singular contributions proportional Ci ej. Charge conjugation 
implies: 

p(i) ^ p(i) ^ p(i) ^ p(i) 
^qjqi^qi ^qjqt^qi ^g^gi^gi ^qjqi^qi 



(26) 

qjqj^qi ^ " qjqj^qi ~ " qjqj^qi ~ " qjqj^qi 



p(i) ^ pW ^ p(i) _ ^ p(i) 



The cross sections (Jq^q^ , with the same quark flavor in both the initial and final state, can also be classified according 
to the charge factor that accompanies the simple pole. For the case in which the charge factor is with j ^ i (that 
is contributions coming from the first two diagrams in Figure ^ , the factorization of forward singularities fixes the 
sum of Pqlq-^-q- (m, v) and Pq^^^^q. (w, w), but can not determine them individually. However, the symmetry under the 
exchange of the quark-antiquark pair of flavor j of the cross section, implies 

p(i) _ p(i) (27) 

^q^qj^qi — ^qtqj^qi y v^w 



whereas charge conjugation gives 



qiqj^qi ^ qiqj*-qi ■ \ ' 



Again there are no singularities in the e^ej pieces. 

The terms proportional to e| in the cr^ig; cross section contribute to the combination of the Pq^<—q and ^'g^Lg 
kernels. In this case, although charge conjugation implies 



(1) _ p(i) 

(1) _ p(T) - (29) 



py-) _ TDK'-) 

^qq^q — ^qq^q 

p(l) ^ p(l) 

qq^q qq^q 



there is no symmetry at the cross section level relating Pqq}—q and -Pg^Lgj ^^nd thus allowing to disentangle them. This 
means that the cancellation of the forward singularities in this particular process would not allow to obtain the scale 
dependence of the valence non-singlet combinations of quark fracture functions. 

On the other hand, singlet combinations of fracture functions, and also those non-singlet combinations that include 
Mq,h/p + Mq f^^p, evolve with precisely the above mentioned combination of kernels. Notice that these last two 
combinations are the ones that actually occur in the hadronic cross sections, thus, although the kernels obtained so 
far are not enough to perform the evolution of individual flavors, they suffice to obtain the scale dependence of the 
cross section. 

Finally, the simple pole in the Cgg cross section is factorized into the Pjjq^^q kernel, which due to charge conjugation 
satisfies: 

3(1) ^ p(i 
qq^q qq 

Explicit expressions for the above mentioned kernels are given in Appendix B 



P^lq = P^l,- (30) 
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Having computed the relevant kernels, the evolution equations for the corresponding fracture function flavor combi- 
nations can be obtained taking moments in two variables of their respective factorization prescriptions, equation l|20|l . 
as it is done in reference jl3| • Taking the derivative with respect to the scale , the resulting algebraic equations for 
the fracture functions moments, can then be solved using the already known explicit solutions for those of the parton 
densities and fragmentation functions. Inverting moments and using eqs. (17-20) wc finally arrive to: 

where, again, we have chosen all the scales to be equal. At order as, the Dirac delta of the non homogeneous kernels, 
equation (|23|l . allows to trivially perform the integration over the v variable for the non homogeneous terms in equation 

(|31|l . which can be recasted in the more familiar form of reference 0|. However, at NLO the kernels Pj.l^j{u,v) are 
no longer constrained to the curve v = {1 — u)/u and thus the general expressions for the evolution equations involve 
a double convolution. 

IV. PHENOMENOLOGY 

In this section we evaluate the phenomenological consequences of the non homogeneous corrections computed so 
far, examining the impact of these corrections relative to the homogeneous evolution pol |. Unfortunately, there is 
scarce direct information on semi-inclusive cross sections with a wide kinematical coverage, as required as input for 
the evolution equations. However there are sensible model estimates available, which can used for the present purpose, 
in particular for tt"*" production from proton targets. 

In order to implement the comparison, we consider a model estimate of the corresponding fracture function M^^+ /p 
based on the ideas of reference |2j| for the description of forward hadrons in DIS, already applied to fracture functions 
in reference 0- In this approach, the fracture function is written as 

Mq^^+/p{xB,XL,Ql),^(l)^+/p{xL)q{xB,Ql) , (32) 

in terms of the flux of positive pions in the proton times the quark densities in a neutron, x^ is the ratio between 
the final state hadron energy and that of the proton target in the laboratory frame. For very forward hadrons, 
XL {1 — XB)vh- Taking this estimate as an adequate approximation for Mq .^+/p at given initial scale Qq, it is 
possible to obtain the fracture densities at higher scales using the evolution equations obtained in the last section, 
and realistic parameterizations for parton densities |22j and fragmentation functions 23] . 

As we have anticipated, although the results obtained in the previous section, together with those of reference 
|lO| for the gluon initiated contributions, complete the QCD corrections up to order-a^ in the cross sections and 
thus the NLO kernels required for Mq, we still do not know those required for the evolution of the gluon fracture 
density Mg. These kernels show up in the order-a^ cross section and contribute indirectly to the evolution of Mq 
through the homogeneous terms in the evolution equations, meanwhile it is worth assessing the effects of its direct 
NLO contributions. 

In Figure |S1 we plot the outcome of evolving the model estimate for Af^ ,r+ /p from an initial scale Qq — 1 GeV^ 
to — 100 Gey^ as a function of a;^ for four different values of x^. The solid lines correspond to the full NLO 
evolution for the flavor singlet combination of Af^ ,r+ yp 

<»iL^E^-^..-VP' (33) 

the dashed lines come from using LO evolution, and the dots show the evolution using only homogeneous terms. 

NLO non homogeneous corrections clearly have a different relative impact depending on the kinematical region 
analysed. For example, in the low xl region, they dominate over the homogeneous terms and their effects become as 
large as those of the LO terms for large xb- In this region the difference between LO and NLO evolution is almost 
independent of a;^, however, the LO non homogeneous evolution effects increase considerably as a;^ falls. In Figure 
Elwe show the same as in Figure |31 but as a function of for = 0.15. 

For larger fractions of hadron momentum, both the LO and NLO non homogeneous contributions become more 
and more suppressed, and although NLO effects overpower LO corrections, they can not match the growth of the 



diogA'P 27r 
as(M^) 1 /-sTc dM 



du 
u 



pi!_l-(.) + ^/^«(.) 



27r 



dv 

V 



5(0) 



Pr 



^sjM') p(l) 



2tt 



ki^j 



{u,v) 
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FIG. 5: Flavor Singlet combination of xMq at — 100 GeV^ as evolved from Qq — 1 GeV^ including NLO non homogeneus 
corrections for different values of the pion energy fraction xl- Solid lines correspond to the full NLO evolution, dashed lines 
correspond to LO evolution and dotted lines show the evolution using only homogeneous terms. 



3 I . — ^ — I 0.6 r ' — ' — I 0.29 r ' — ' — I 0.2 




FIG. 6: Scale dependence of xMq for xb = 0.15 



homogeneous contributions. Large hadron momentum fraction is precisely the region relevant for hard diffractive 
interactions. The lack of significant non homogeneus effects in the evolution in this region, is consistent with previous 
studies on the scale dependence of diffractive deep inelastic data . 

Notice that the relative importance of homogeneous and non homogeneous contributions to the evolution depends 
crucially on the the size and shape of the appropriate fragmentation and fracture functions corresponding to the process 
considered. This means that the specific energy fraction at which non homogeneus effects become non negligible is 
characteristic of each process, and in this particular case, depend on the model estimate used for the fracture function, 
as well as on the parameterizations used for parton densities and fragmentation functions. From another side, the 
value chosen in the present estimate as departing point for the evolution, i.e. the value Qo *t which we assume the 
model estimate to hold, modifies to some extent the outcome of the evolution and thus the onset of non homogeneous 
effects. In any case, the numerical estimates shown in this section illustrate the sort of scale dependence one can 
expect from a given one particle inclusive processes in the light of our present knowledge of QCD corrections. 

V. CONCLUSIONS 

We have computed the 0{a1) quark initiated corrections to one particle inclusive deep inelastic scattering cross 
sections, verifying and extending for these processes the feasibility of the approach proposed in reference ^3 to deal 
with overlapping singularities. Although quark initiated corrections include a more complex singularity structure, 
once backward configurations are regulated, forward singularities can be treated in a similar way to those in gluon 
initiated corrections. 

We made an explicit check of the factorization of coUinear singularities arising from quark initiated corrections 
at this order and obtained the evolution kernels required to compute non homogeneous contributions to the scale 
dependence. These kernels, which are labeled by three indices, corresponding to the nature of the parton that 
interacts with photon, hadronizes, or initiates the hard processes, respectively, show up entangled in the cross section 
and can not be extracted individually for each flavor from SIDIS. However, flavor combinations required for the singlet 
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and total flavor combinations can be obtained. Finally, using a model estimate for fracture functions of protons into 
pions, we found that the impact in the scale dependence of non homogeneous NLO corrections relative to LO and 
to homogeneous terms depends strongly on the kinematical region under analysis. Specifically, NLO corrections are 
particularly large for small values of parton and pion momentum fractions. 
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Appendix A 



In this appendix we present the angular integrals required in the computation of the quark initiated corrections to 
SIDIS at order 0{al). As usual, all these integrals can be put, by means of partial fractioning, in the standard form 

M 



/(fc,0= / dPi 
Jo 



sin^+' (3i sin' (32 



df32 



(a + 5cos^i)'=(A + Bcosf3i + Csin/3i cos/32)' ' 



(34) 



and can be classified according to whether their parameters satisfy either — or = + C^, both relations 
simultaneously, or neither of them. In the following table, these cases are labeled as Ia{k,l), lA{k,l), lAa{k,l) and 
lQ{k,l) respectively. The parameters a, 6, A, B and C, are functions of the variables x^, u, v and w. The resulting 
integrals as well as their coefficients can develop further singularities associated to these variables, what requires some 
care when expanding in power series of the regulator e. Specifically, the limit w —> 1 implies always that b — > ±a, 
C ^ and B — > ±A, which means that factors like {Ah — a BY or (a^ — 6^)' often present in the expressions have to 
be kept to all orders. This factorization is trivial for the integrals that can be computed to all orders in e, but as it was 
mentioned in Section II, some special care is required for the integrals that have to be expanded before calculation. 
All the integrals of this last type can be handled with the method described in reference pX]| . 
The integrals that can be done to all orders in e are: 



• Iq: ^ 6^ and A^ ^ B^ 

2n 



C2 



lo 
lo 

lo 
lo 
lo 
lo 
lo 
lo 
lo 
lo 
lo 



0,0) 



1 + e' 
1,0) = a/o(0,0) 



-1,-1) = 
-1,-2) = 
-2,0) = 

-2,-1) = 
-2,-2) = 
-3,0) = 
-3,-1) = 



a A 4 
aA^ 



bB 



/o(0,0). 



3 + e_ 
2ABb + a{B^ + C^) 



a' A 



Ab-" 



3 + e 
/o(0,0), 
2abB 



4(0,0) 



, 3a52 



/o(0,0), 



AH^ + AAaBb + a^{B^ + C^) b^3 B^ + C^) 



4(0,0), 



3ab{Ab + aB) 



Zb^B 



(3 + e)(5 + e) 



27r 



l^Q)^ — \l-{a + by 



e 

2' ' 



(3 + e)(5 + e) 



4(0,0), 
2b 



4(0,0), 



1,-1) 



Ab-aB 



4(1,0) 



2ttB 



(35) 
(36) 
(37) 

(38) 

(39) 

(40) 

(41) 

(42) 

(43) 

(44) 

(45) 
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/o(l,-2) = 
/o(l,-3) = 
+ 

/o(l,-4) = 



1 
1 

63 



{Ab-aBy - 



2 C2(a2-62M 



2 + e 



/o(l,0) 



2tt 



2ABb-aB^ 



aC 



2 1 



2 + e 



{Ab-aBy - 
2n 



3 3{Ab-aB){a^ -b^)C^ 



2 + e 



/o(l,0) 



(1 + e) 63 
1 



. -w., X 9 9s b'^B^ ^9 f Bb"^ a(Ab-aB) 

B{3AbiAb-aB)+a^B^) + ^^—+3C^' ' ^ ' 



3 + e 



64 



{Ab-aB) 

271 



3 + e 

2\2 ri2n 



2 + e 



(46) 



(47) 



2 + e 

4ABb{Ab-aBf + aB^{2A'^b^ -a^B^ 



/o(l,0) 



(2 + e)(4 + e) 

{iABb- aB^)b^ B 



(1 + e) 64 

6C^b'^B{Ab~aB) &C^a{Ab-aBy 3aC^{a^-b^ 



3 + e 
3a6C2 



/o(2,0)- 

/o(2,-l) 
/o(2,-2) 

/o(2,-3) 



a e 



a2 -62 
1 

eab 
1 



3 + e 



/o(l,0) 



2 + e 



2tt 



i^~b^ 



[a{Ab- aB)e - B{a^ - b^)] Io{2, 0) + 



(2 + e)(3 + e) (3 + e)(4 + e) 



2nB 



27r 
eab"^ 

1 
e a 



63 



[a(A6-aB)2e-2B(a2 -62)(^6-aS) -aC2(a2 -62)] /o(2,0) 
'B{2Ab-aB)+ ''^^l-p ' 

a{Ab - aBfe ~ -iB {a^ - b^) {Ab - aBf - ■iC'^ {a^ ~ b 



(48) 

(49) 
(50) 

(51) 



X \a{Ab-aB) 



B (a2 - 62 
2 + e 



/o(2,0) 



27r 

£ a 63 



3C^ 



/62b a(A6-aB) 



a2B 



(l + e)(2 + e) 



B ['iAb{Ab-aB)+a^ B"^ 



2a2 B2 -3a^5B\ 



1 + e 



■J 



(52) 



where a > 6 > 0. Notice that the cases — b^ and A^ ~ B^ + corresponding to integrals (|36|l - 143() can be 
easily obtained taking the appropriate limit, whereas for integrals (|44|I - H52|I only the limit 6*2 ^2 — B^ can 
be taken, giving the corresponding I a in tegral. The hypergeometric function that appears in the cases fc > 
can be expanded in power series in e |24| : 



2F1 



-,e,l + e;z 



2 3 

= 1 + 1- Li^iz) + ^( log2(l - z) log(z) + 2 log(l - z) Li2(l - z) 
- 2 Lis (1 - ^) + 2 C(3) - 4 Li3(z)) + 0(e4) 



(53) 



lAa'- a? — b^ , A^ — B^ — = 0. Using the results in |2^, all these integrals can be put in the standard form 
I2I: 



lAa{k,l) 



2n 



a '^-^3^-^±B(l-k+i,l-l + i 



r(i + f)2 



2F1 



k.Ll 



e A-B 



2' 2A 



(54) 



where < — ^, -j < 1. The hypergeometric hmctions appearing in the cases k,l > should be carefully treated 
as they usually develop singularities in the limit B — > —A. Using the relations in sections (2.8) and (2.9) in 
the above mentioned hypergeometric functions can be written always in terms of 



2J^l 



e e ^ e A-B 
2'2' ^ 2' 2A 



(55) 



which is regular when B = — A and admits the expansion [2 

2-Fl 



111 ' 
2' 2' 



e2 e3 

1 + - Li2{z) + Yg ( log'(l - z) \og{z) + 2 log(l - z) Li2(l - z) 
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2 Li3(l - z) + 2 C(3) - 2 Lis (z)) + 0{e^) 



(56) 



All the integrals which cannot be calculated to all orders in e are of the classes /^(fc, or /a(fc, with k,l > 0. With 
an adequate choice of the reference frame, the latter integrals can be recasted in terms of lA{k, I) integrals also, which 
are considerably simpler to calculate. The 10 resulting integrals belong to the sets Ia{^,)\ ^^(2,1), /a(1,2) and 
/^(2, 2). In each case, the potentially singular factors have to be extracted, as described in before the expansion 
in powers of e. With the exception of two integrals of the type Ia{^, 1), all the integrals present potential singularities 
due to a factor {Aa — with n = 1,2,3. In the remaining two, corresponding to the 7^(1,1) case, a factor 

(a^ — b'^y/^ has to be kept to all orders. Explicitly: 



/a(1,1) = Tr{Aa~bB) 

a + b 



log 



log 



6r(i + i) 

b{A-B) 
A{a-b) 



1 1 1 +1 HA~B) 
2'2' 2' A{a-b) 



log 



A (a + b) 
Aa-bB 



log 



f A{a + b) 
[Aa-bB 



log 



A (a + b) 



62 {A- B)^ (Aa-bB) 



47r 



bjA + B) 
Aa-bB 



Li2 



/Aa-bB 
\A{a-b) 



+ 0{e^) ) , case 1 



t 

X {a-b)-''^ [Aa-bB)-'''^2Fi 
log^ {A [a - 6)) - 2 log 



2' 2' 



e b{A + B) 

2'"' ' Aa-6B 



A{a-b) 



7r(a2 -62)^/2 [2 



log 



^ (g - 5) 



A 



log 



A (g - 6) 
^g - 



log^ {Aa-bB) 



0(e2) I , case 2 



/a(2, 1) = (g + 5)^1+^/2 |(A + B) (Aa -bB)-'-'/^ 2F1 



e fe(A + S) 

2'''^ + ''"Ag-6S 



2(a-6)-i-^/^Fi 



e ^ 26 
2 g — 



2tt 



TT r(i 



X (Aa-6S)-2+^/2 (2-e)2Fi 

-^(Aa-6B)-2+^/2 
6 (A + S) 



2' 2' 2'"g 



{a -b) (Aa-bB) 
26 



e rfi 



£66 &(A-i?) 

2'2' 2' A(a-6) 



/ 6(A + g) 
I g + 6 



A log 



/ A (g + 6) 
V Ag-6B 



A{a-b)iA + B) 
(a + b){A- B) 



log(g - 6) 



, (A + B){Aa-bB) ^ . 
log(yl) - — ^ log(g + 6) + A log 



g + 6 ' {a + b){A-B) 

2 {A{Ab - a B) + {A - B) {Aa - b B)) 



a — b 



log - 



b{A-B) 
A{a-b) 



2 log 



Li2 - 



(g + 6) (A - B) 

A (a + b)' 



log 



- 65 

26 



g 



log 
0(e2 



A (a + 6) 



62 (A- B)^ (Aa-bB) 



A (a + b) 
Aa -bB 

A( Li2 



Aa -bB 
A{a-b) 



2U2 - 



b{A + B) 
Aa-bB 



Ia{1,2) = - 



r(i + e) 



(Aa-6B)2|4r(l + i)2 

6 (a -6) (A + B) (2-e) (4 - e) 
Aa -bB 



(g + 6)(2 - e) + ^ (a - b)-'/\Aa - bBf'^ (^2 (g + 6) e - 86 

e e b{A-By 



2F1 



A{a-b) 



1 
A 



2{Aa b B) _ 1 5 + 5)^./2 + 5)./2(^ ^ _ ^5)-e/2 



2 6B- A(a-6) 
b{A + B) 



2 Aa-bB 



(57) 



(58) 



(59) 
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+ 



4:A{Aa-bB) 

A{a-b) + 2Bb /A{a + b) 



AA 
b{Ab-aB) 



2A(Aa~bB) Aa-bB 
log{A{a~b)) 

log 



biAb-aB)^^j Aia + bl ^^^ 



Aa -bB 



Aa + 2Ab + Bb ^ fa + b 
— log 



A{a^ + b'^)-2abB 



4:{Aa-bB) 



4A 

b{Ab-aB) 



A{Aa~bB) 
' Aa- bB 



\Aa~bB 



log 



a — b 

a + b 



b(Ab-aB), /A (a + b) 
log 



+ Li2 



A{a-b) 



Aa -bB 
hO(e2) 



2{Aa-bB) ' ° Va-6 
A (a + b) 



log 



Aa-bB °\Aa-bB 
b{A-B) 



A{a-b) 
b{Ab-aB) 

b^ {A- B)^ (Aa-bB) J ' 2{Aa-bB) 



biA + B) \ 
Aa-bB 



(60) 



/a (2, 2) 



r(i 



^)(" + ^) {Aia-b){4-e)-2{Aa-bB)) + - A~'/^ {a - b)''l^ 
a — b e 



(Aa-6B)3 [gr(l + l) 

< (A a - 6 ( 4 ^2 &2 (3 _ g) (4 _ g) _ 2 e (2 - e) (A a + 6 + A 6 (A (a - &) - S (a + 6)) e (4 - e)^ 



+ 2AaBb{2-e)e'^ -ABb{l-e){'i-e){Aa{2-e) + Bb)]2Fi 
4-e 



1 1 1+1 M^-g) 
2'2' 2' ^(a-6) 



4(A(a + 6)-3&(A + 5))-^i4^^^ + 2KA + 5)e. ^ + 



2-e ■ " ' ' ~ ' A{a + b) 
X ^ 5)<:/2 a _ ^ 5)-^/2 [2A{a + b){2-() + b{A + B) 6)2^1 



1 + 
2'"' ' Aa-6B 



(Aa-6B) 



= /2 



X loa 



(A{a 


^b)\ 






\Aa- 
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+ e 
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2(a + 6)(Aa-6S) ^ - 65 

2 62(yl-S) 2 62(A + B)2 2 62(^2-^2 



A{a-b) 



a — b 



b^ {A + B) (5 {Ab-aB) + 3{Aa-bB) 



A {a + b) {Aa-bB) 



log {A {a -b)) 



A{a + b) Aa-bB 

6b{Ab- aB) - {a + 3b){Aa-bB) 



X log 
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b{Aa-bB + 2{Ab-aB) 
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log 



A{Aa-bB) 
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Aa- bB 



log 



A{a-b) 
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X log ( + b{?.A{Ab-aB) + B{Aa-bB)) / A{a + b) 



A{a-b) 



A{Aa-bB) 



b{SA{Ab-aB)+B{Aa-bB)) f ( b{A + B) 

L12 



2{Aa-bB) 



Aa -bB 



Li, 



Aa- bB 

Aa-bB 

A{a-b) 



a — b ^ 
A {a + b) 



p (A-B)2 {Aa-bB) 
f 0(e2 



(61) 
(62) 



For the Ia{^, 1) integrals, case 1 means that in the limit w —^ 1 the parameters satisfy b — > —a and B —>■ —A, whereas 
in case 2 they satisfy 6 —a and B ^ A. In this last case the factors {Aa — bB)~"' remain finite when w ^ 1. 

Notice that in the above list, there are some integrals that can be obtained from another by a suitable transformation, 
for example, integral /o(2, 0) can be obtained from /o(l, 0) by derivation with respect to parameter a. Further relations 
can be obtained using the results presented in [23 |. 



Appendix B 



In this appendix we present results for the NLO non homogeneous kernels for fracture functions obtained from the 
factorization of coUinear singularities in the 0{a].) SIDIS cross section. 
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In the case where a gluon is the fragmenting particle, we obtained 
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1 + X 



Pqg^q{x) = y— ^ , and a/ = 



l-u 



(63) 



X u 

For the kernels relating different quark species, comig from the quark fragmenting processes, we obtained 

[l+u + u^] v-u^ [l-v + v"^) - (l -u + 'u^) [l+Auv + 2uv'^) +2pg^g{u) 
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P^jl^g^{u,v) = CFTF{2u 
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u{l + vY 



(65) 



+ ( log(u) + log(a/ -v) - \og{v) 
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+ 2 log(l + v)) 



u{l + vr J 
u + v'^ + 7uv {1 - u - uv) - {1 - u) (l + v + uv + u^v^) 



+ 



{1 — uv) 

1 / \ 1 / x 1 / x ^ f 1 — U — UV 

log(u) + log(a/ -v)+ \og{v) - 2 log I 



M (l + U^) 

(1 — ut;)^ 



(1 + V - 2 (1 - u) u (1 + v)) ) , 



(66) 



(67) 



where 

Pq^g{x) = l-2x + 2x'^. 
Finally, for the kernels relating quarks and antiquarks of the same flavor, 



(1) J 2u(l + 2u2) (4 + 3^^ + 5^2) w 2u^ {3 + 2u + 3u'^) v 

Pqq^q(u,v) =CfTf { — -^2 ^ 



3(l-u)' 



3(l-u)' 



8(l + u) 2 (3 + 10u-3u=^) 4(l-2u) 

"I T ! T To I r. I-, : ^ I" 



u{l + vf Zu{l + vf 



3 (l + v) 



{1-uY u{l + vy 

2u (5-3u + 4u2) 



3(l-u)^ 

log(l - u) 

2u(5-6u + 4u2) Au^ {2-u + 2u'^) v 2u^ {l + u^) v'^ 4 



2u'^ [7 -2u + 7u'^) V 4u3(l + u2)t;2 2(l + «2) 



3(l-< 



{l-uf Z{l-u){l + v) 



3(l-u)^ 3(l-u)'' 
4(l + u) 2(l + u)^ 4(2-7/2) 



+ 



u{l + vf u{l + vy 3(l-u)(l + u) 
4^2 (2-M + 2u2) 2u^ {l+u"^) v"^ 



3(l-u)' 
4 (2 - 



+ 



log(u) + 
4 



{l-uf u{l+vy 
2u (5 - 6m + 4u2) 



3(l-w)' 
4(1 + 1*) 2(1 + 1*)^ 



3 (1-u) (l + v) 



(l-u) 
log(a/ - v) 



u{l+v) u{l + v) u{l + v) 

2u2 2u2 (l +^2) ^ ^ (l +m2) i;2 



{i-uY {i-uY 



+ 



4(l + u) 2(1 + m)^ 2(1 + m) 



+ 



u{i + vf u{i + vy u{i + vy i + v 

\2 



log(t;) + 



(68) 



(1 - uf 
2(1-2u)m_^ 2u'^v 



3 (1 -m) 3 (1 -«) 



+ 



{1 + u) ^ 4(l+u)^ 2 (5-7u2) 



^.{1 + vf u{l + vf u{l + vf 3(l-u)(l + t;) 



log(l + v) 



(69) 



and 



P^llq{u,v)+P!,^qlq{u,v) 



CfTfI - 



2 u (2 + u2) 4 u2 (5 + 3 u + 4 u2) 2 (3 + 2 u + 3 ^2) 



3 (l-u)' 



3 (l-u)' 



{1-uy 



u{i + vy 



8 (l + u) _ 2 (3 + 14u-3u2) ^ 4 (_2 + u) ^ IQu^ 16«2(i + u) 
^~ ~ ~~ ~' 3(1 + 1;) ~ {1-uvf^ {l-uvf 



u{l + vf 2,u{l + vf 

4u(l + 4u-u2) 2^ (5- 12w + 67i2 



{1 — uv) 



+ 



3 (1 - u) (1 - uv) 

2 u2 (5 + 2 M + 5 ^2) ^; 4 (1 + m2) v'^ 



2 u (4 - w + 3 1 

3(1- w,)' 
2u (1 + m2) 



3(l-«)' 



(1-u)^ 3 (l-«) {l-uv) 



log(l - u) 
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+ 



2u{2 + u) 2u^ (l + Au + u^) V 2u^{l + u'^) 



3{l-uy 



3(1 -u) 



+ 



(1 - uY 



4 4(1 + u) 

U (1 + w)^ M (1 + vf 



+ 



+ 



2{l + uf 2{l+u) 



(i + u) 4m(1 + m)^ 



u{l + vf 1 + v 

2u (5-8w2) 
3 (1 -u) {l-uv) 



{1 — uv) {1 — uv) 



(1 — u vY 



log(u) + 



2u' 



{l-uY 



2v? {l + u^) V ^2 (l + i;2 



(i-«r 



4(1 + ;/) 2(1 + (()" 



u{l + vY u{l + vY u{l + vY 



4u(l + u) 



(! + ") 
1 + v 



+ 



2 ?i (5 - 6 1 



(l-?i«)2 3 [l-u) [l-uv) 



log(a/ - - 



4 (i + u + ^; 2 (l + V?-) V 



3{l-uY 
2{l+u) 



+ 



(1 - uf 

8 u- (1 + h) 
(1 — ut;) [l — uvY 
2 u (4 - 4 u + 3 u^) 

3(1-m)' 

4(l + u) 2(l + u)^ 



+ 



u{l-\-vY u{l + vf u{l + vY 



{i-uY 

iu^il + u) 4m(1 + m)^ 2m (8-5^2) 



l + v 



{l-uvY {l-uvf {l-uvY 'i {l-u) {l-uv) 



\og{v) 



+ 



2u + 



8 8(1 + m) 4(1 + ^)^^ 2(7+12m + 7m^) 

« (1 + t;)^ ~ u (1 + vf u{l+vf 3 (1 + u) (1 + v) 



2u{l + v? 



+ 



3 (1 + u) [l-uv] 

16^3 1«..2 



log(l +v) + 



2{h-Qu)u 2u'^v 



3 (1 -u) 

2 



[l — uv) [l — uv) 

28 u (l + 2?i + u2) 
3 (1+w) (l-uv) 



16u2(l + u) 8u(l + u) 

"I 7" To I" 



2 (1 + M^) 

3 (l-u) 3(l + u)(l + t;) 

4n2 



(1 — uv) 

log(l — uv) > . 



3 (1 + m) {l-uv) 



(70) 



As mentioned in Section III, in this last case, only the combination Pgg]_q(w, v) + Pql}—q{u, v) can be extracted. The 
remaining quark initiated kernels can be obtained from the charge conjugation relations given in that Section. 
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